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Abstract 

Selection on the level of loosely associated groups has been 
suggested as a route towards the evolution of cooperation be- 
tween individuals and the subsequent formation of higher- 
level biological entities. Such group selection explanations 
remain problematic, however, due to the narrow range of pa- 
rameters under which they can overturn within-group selec- 
tion that favours selfish behaviour. In principle, individual 
selection could act on such parameters so as to strengthen the 
force of between-group selection and hence increase coopera- 
tion and individual fitness, as illustrated in our previous work. 
However, such a process cannot operate in parameter regions 
where group selection effects are totally absent, since there 
would be no selective gradient to follow. One key parameter, 
which when increased often rapidly causes group selection 
effects to tend to zero, is initial group size, for when groups 
are formed randomly then even moderately sized groups lack 
significant variance in their composition. However, the con- 
sequent restriction of any group selection effect to small sized 
groups is derived from models that assume selfish types will 
competitively exclude their more cooperative counterparts at 
within-group equilibrium. In such cases, diversity in the mi- 
grant pool can tend to zero and accordingly variance in group 
composition cannot be generated. In contrast, we show that if 
within-group dynamics lead to a stable coexistence of selfish 
and cooperative types, then the range of group sizes showing 
some effect of group selection is much larger. 



Introduction 

The evolution of cooperation between biological individ- 
uals, both generally and as a vital part of the formation 
of new higher-level composite individuals, is an important 
and m uch discussed open question in both evolutionary bi- 
ology dMavnard Smith and Szathmarvi 1 19951: iKelleri [l999 ; 



Michodl 1 19991: iHammersteini [20031: lOkashai 120061) and ar- 
tificial life (BedauetaUioOO). The fundamental problem 
is that any group of cooperative types, whose members do- 
nate some component of individual fitness in order to ben- 
efit their group, is vulnerable to invasion by selfish cheats 
that reap the group benefits of cooperation without paying 
the individual cost. Various mechanisms by which cooper- 
ation can nevertheless evolve have been suggested, includ- 
ing so-called 'group selection' models in which population 



structure exists such that individuals spend part of their time 
in groups (rather than freely mixed in the whole population) 
before mixing in a m igrant pool from which new groups are 
formed dWilsonl ll9801. Although selection within any given 
group will always favour selfish individuals, groups with a 
higher proportion of cooperators are more productive and 
hence contribute more individuals to the migrant pool and 
the next generation of groups. 

Such models have been found to allow cooperation to 
evolve or be preserved in a population, but only within 
certain narrow parameter ranges (discussed below). This 
clearly presents a problem when appealing to such mech- 
anisms as a route towards the evolution of hig her-level indi- 



vidua ls. As we have argued in previous work ([Powers et al 



2007), the conditions that allow group selection to be effec- 
tive and control its strength need not be externally imposed. 
Instead, they may be products of individual characters and 
hence subject to individual selection. Specifically, if key pa- 
rameters such as group size are subject to individual adapta- 
tion (for example, via production of extracellular matrix in 
a bac terial biofilm), then a pr ocess akin to niche construc- 
tion (Odli ng-Smee et al. L l2003h supporting the evolution of 
cooperation may occur whereby cooperative traits and those 
affecting the strength of group selection evolve concurrently. 

However, this process would require the existence of an 
adaptive gradient, such that a small parameter change could 
increase the strength of group selection and consequently 
cooperation and individual fitness. In this paper we con- 
sider one important parameter, group size. Where groups 
are formed randomly, existing models have shown that in- 
creasing group size rapidly c auses the me asurable effect of 
group selection to reach zero dWadelll978h . This then means 
that for a large region of parameter space, a small decrease in 
group size, e.g. by individual mutation, would have no effect 
on fitness, i.e. there would be no selective gradient towards 
smaller groups and increased individual fitness. However, 
the model developed and presented in this paper suggests 
that the rapid tendency of group selection effects to zero is 
a consequence of an assumption of directional within-group 
selection. Specifically, classical group selection models at- 



tempt to explain the global promotion of a cooperative allele 
that is driven extinct at within-group equilibrium through 
directional selection for a rival selfish allele. While this 
assumption of directional selection and hence competitive 
exclusion of types is commonplace in population genetics 
models, a competitive coexistence of two types is instead 
often permitted in ecological models. For example, the clas- 
sical Lotka-Volterra competition equati ons allow for coexis- 
tence as well as exclusion dMavUl976l) . Such dynamics are 
of relevance to potential multi-species group selection sce- 
narios; for example, during egalitarian major evolutionary 
transitions in which different unrelated individuals eventu- 
ally form a new level of selection, and with in extant multi- 
speci es consortia such as bacterial biofilms dBurmolle et all 
2006) in which group selection effects may be pertinent, [j 
Results presented in this paper show that where a within- 
group stable coexistence of types exists, measurable group 
selection effects are sustained over a much larger range of 
group sizes, potentially providing an individual adaptive 
gradient towards smaller groups that enhance group selec- 
tion over a much larger range of parameter space. 

The Limits of Group Selection 

Group selection theory is typically understood as the idea 
that the differential productivity of groups can lead to 
changes in allele frequency in a gene pool, in a manner anal- 
ogous to individual selection acting through the d i fferen - 
tial productivity of individuals dWadd, 1 1 978t IWilson , 1980). 
In this paper we are concerned with type 1 group selec- 
tion, which defines gr oup fitness as the average fitness of 
the group members^ (O kasha , 2006, chap. 2). Models of 
this process therefore seek to investigate the effect of group 
structure on the evolution of individual social traits such as 
cooperation. In particular, a trait that is individually dis- 
advantageous may nevertheless evolve if it has a positive 
effect on the g roup as a whole, an idea that dates back to 
Darwinl dl87lb . Using the language of multilevel selection 
theory, such cooperative traits are selected against within- 
group, since they confer a fitness disadvantage relative to 
other group members, but are favoured under between-group 
selection, since they increase group productivity, i.e. aver- 
age absolute group member fitness. Whether or not the trait 
spreads in the global gene pool then depends on the balance 
of these two selective forces. 

Three factors are pertinent in determining the outcome 
of such a scenario. The first is the individual cost to 
group benefit ratio of the cooperative act; a large individ- 
ual cost strengthens within-group selection against coopera- 
tion, while a large group benefit strengthens between-group 



selection. Secondly, a group mixing mechanism must ex- 
ist so that the increased productivity of more cooperative 
groups may affect the global allele frequencies. We consider 
here a multi-generational variant of D.S. Wilson 's dl980h 
trait-group model, where a global mixing stage occurs every 
t generations in which the progeny of all groups disperse, 
join a global migrant pool, and then reform new groups of 
random composition. Since cooperative groups will con- 
tain more individuals prior to dispersal, they will constitute 
a larger fraction of the global migrant pool, thereby bias- 
ing the global allele frequencies. How often global mix- 
ing occurs is therefore a crucial parameter, since groups 
must be mixed before the within-group equilibrium has been 
reached, otherwise there will be no difference in group pro- 
ductivity for selection to act on (assuming selectio n within 
all gr oups leads to the same group equilibrium d Wilson , 



1992)) 



'Coexistence dynamics may also apply to single-species sce- 
narios under certain regimes such as balancing selection (see be- 
low) and are hence relevant to group selection discussions more 
generally. 

2 Type 2 defines group fitness as number of offspring groups. 



The third factor, and the one which is most often com- 
mented on in the group selection literature, is the require- 
ment for there to be variation in the groups' allelic compo- 
sition. The larger this variance, the greater the difference 
in group productivity (since within-group dynamics are as- 
sumed to be deterministic) and likewise the effect of group 
selection. Many theoretical models have shown that if group 
composition is random then very small initial group sizes are 
needed in order to produce the betw een-group va riance nec- 
essary for a measurable effect (see dWadeUl978l) for a clas- 
sical review). It is therefore usually concluded that some 
kind of non-random group formation is required in order 
for group selection to have any significant effect. The most 
common way that such assortative grouping is believed to 
occur in nature is through kin grouping, where the group 
members are related by descent from a common ancestor 
and are hence more similar to each other than t o members 
of oth er groups. Consequently, kin selection dHamilton , 
Il964allbl) is commonly seen as the only pertinent force in 
social evolution. However, even with only simple random 
sampling, the range of group sizes which produce a group 
selection effect may be strongly affected if the assumption 
of within-group competitive exclusion is relaxed. This is 
because possible between-group variance from sampling er- 
ror is controlled by the frequency of the least frequent type 
in the migrant pool. In the competitive exclusion case the 
possible variance rapidly tends to zero as the selfish type ap- 
proaches fixation. If, however, coexistence of the two types 
is possible at equilibrium then we would expect the possible 
variance to be maintained above zero even at within-group 
equilibrium frequencies. 

Group Competition: Aggregation and Dispersal 
Through a Migrant Pool 

In this paper, we consider a model of a population structure 
where individuals reproduce in groups for a number of gen- 
erations. After this period of reproduction within groups, 



all groups disperse and their progeny mix freely together in 
a migrant pool. Thereafter, new groups are created from 
the individuals in the migrant pool, and the process repeats. 
This proce ss is a multi-generat ional variant of the trait-group 
model of IWilson ( 198dl 1 1987b. and correspond s closely to 
the "Haystack" model of Ma vnard Smifhl d 1 9641) . Examples 
of natural populations that fit this model particularly well 
include soil-dwelling populations of micro-organis ms that 
are oc casionally mixed together during rainstorms dWilson , 
1980, p. 22), and a type of desert leaf cutting ant which lives 




in colonies that periodically disband and whi ch are founded 



by un related females from a mating swarm dRissing et al 
Il989h . However, the issues explored by the general form of 
the model are broadly applicable to a wide range of popula- 
tion structures where individuals have the majority of their 
interactions with a subset of the population. An algorithmic 
description of such an aggregation and dispersal process is 
as follows: 

1 . Initialisation: Initialise the migrant pool with N individ- 
uals; 

2. Group formation (aggregation): Assign individuals in 
the migrant pool randomly to groups of size S; 

3. Reproduction: Perform reproduction within groups for t 
time-steps, as described in the next section; 

4. Migrant pool formation (dispersal): Return the progeny 
of each group to the migrant pool; 

5. Iteration: Repeat from step 2 onwards for a number of 
generations, T. 

Group competition occurs in a population structured in 
this fashion during the dispersal stage, since groups that 
have grown to a larger size will contribute more individu- 
als to the migrant pool. Crucially, this means that coop- 
erative traits that are individually disadvantageous but that 
benefit the group can potentially increase in frequency (see 
FigureQ]), depending on the balance of within- and between- 
group selection. 

A key factor in determining the balance between the levels 
of selection in any mo del of group s election is the variance 
in group composition dWilsonlll980l) . for if there is no vari- 
ance then there is nothing for selection to act on ([Darwin, 
Il859h . In particular, there must be variance in initial group 
composition which causes a variance in group size prior to 
dispersal. In the above model, this variance is generated 
through random sampling of individuals from the migrant 
pool, where the sample size corresponds to the initial group 
size. Since an increased sample size causes a decrease in 
between-sample variance, increasing the initial group size 
decreases betwe en-group var iance and hence the efficacy of 
group selection dWadelll978l) . It is therefore often assumed 
that the upper limits of group size which produce a non-zero 
group selection effect are very small. 
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Figure 1: Cooperation can increase in frequency in the 
migrant pool due to differential group contributions, even 
though it decreases in frequency within each group. 

However, in this paper we present results which suggest 
that this follows from an assumption of within-group dy- 
namics that lead to the competitive exclusion of a coopera- 
tive type by its selfish counterpart at within-group equilib- 
rium. Specifically, we are able to show that where within- 
group dynamics instead lead to a stable coexistence of types, 
then the range of initial group sizes over which an effect of 
group selection can be seen is much larger. This is due to the 
fact that since the cooperative type cannot be driven extinct, 
variance in group composition when sampling from the mi- 
grant pool is always possible. In the next section, we de- 
scribe how both competitive exclusion and coexistence dy- 
namics can be modelled within groups. 

Within-group Dynamics: Competitive Exclusion 
Versus Coexistence 

Classical models of group selection consider a scenario 
where a selfish type ultimately drives its cooperative coun- 
terpart extinct at within-group equilibrium. In partic- 
ular, fitnes s functi o ns of the following form, first pro- 
posed by Wright d 1 9451) but subsequently used in a 
plethora of other models dWill iams and Williams, 1957; 
Mavnard Smithl Il964t ICharnov and Krebsl 1 19751; I Wilson . 
198(1 Il987h . are typically used to model within-group se- 
lection: 



fs = 1 + Pc9 

f c = (l+Pcfl)(l-0) 



(1) 

(2) 



In the above equations, f s and f c denote the per capita fit- 
ness of selfish and cooperative individuals within a group, 
respectively. Cooperators, whose proportion within the 
group is denoted by p c , confer a fitness benefit g on every 
group member. Crucially, both types receive this benefit, 
while only cooperators pay a cost, represented by the selec- 
tion coefficient against cooperation, a. It is then clear that if 



these equations are iterated until equilibrium is reached then 
the selfish type will be driven to fixation within the group. 

Both competitive exclusion and stable coexistence within- 
group dynamics can instead be modelled using the standard 
two-specie s symmetric Lotka-Volterra competition equa- 
tions (e.g. dMavUl976l) ). For implementation purposes, we 
use the following difference equation as a discrete approxi- 
mation: 
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(3) 



In the above equation, N^m is the biomass of species i 
at time t, and Mj the intrinsic per capita growth rate. Each 
species has an intrinsic carrying capacity Ki, which is then 
modified through interspecific interactions. Specifically, the 
per capita effect of species j on species i is given by a^ , the 
coefficient of interaction. All such interactions are compet- 
itive in the above equation, as ensured by the negative sign 
and the stipulation that all a > 0. Similarly, an denotes 
the negative density-dependant effect of species i on itself 
that prevents unbounded exponential growth. This coeffi- 
cient can be seen as representing crowding and can vary for 
different species. 

We define selfish (s) and cooperative (c) strategies in the 
above equation through settings of the within- and between- 
type interaction coefficients. Specifically, a selfish type is 
defined as having a large negative per capita effect on both 
itself (a ss ) and the other type (a cs ). A cooperative type is 
then defined as having correspondingly smaller per capita 
negative effects (a cc and a sc ). A pure group of coopera- 
tors will therefore grow to a larger size than a pure group 
of selfish individuals, creating a group productivity differen- 
tial on which selection can potentially act. However, within 
mixed-groups the selfish type will reach the larger frequency 
(provided that a ss is not too large), since a cs > a sc . In 
other words, cooperators are favoured by between-group se- 
lection, while selfish individuals are favoured under within- 
group selection, exactly as in a classic group selection sce- 
nario. It should be noted that our definition of cooperative 
behaviour corr esponds to weak, rather than strong, altruism 
dWilsonl H9801 This follows because although cooperation 
confers a relative fitness disadvantage compared to a selfish 
individual within the same group, it nevertheless increases 
the absolute fitness of all group members, including the co- 
operator. 

It is well known that a stable coexistence of both types 
occurs in such a model when competition for resources 
(space, food, etc.) between individuals of the same type 
is stronger than competition between individuals of differ- 
ent types. Such a case corresponds at the ecological level to 
species occupying different niches, i.e. only partial ly over- 
lapping in their resource requirements (.May., 19761) . Con- 
versely, if between-type competition is stronger than within- 



type competition then competitive exclusion of one type will 
occur. Between- and within-type competition are both mod- 
elled in the Lotka-Volterra equations through the settings of 
the interaction coefficients. Throughout this paper, we as- 
sume the following: 

1. a cc < a ss and a sc < a cs , i.e. that cooperators have 
lower negative density-dependant effects on themselves 
and others; 

2. a cs a sc < 1; 

3. M and K, the intrinsic per capita growth rates and carry- 
ing capacities respectively, are the same for both types. 

Given these assumptions, competitive exclusion of the co- 
operative type occurs when a ss < a cs , producing qualita- 
tively similar dynamics to those of the traditional within- 
group selection equations (Q~|) and |0. However, when 
a ss > a cs then the cooperative type is maintained at within- 
group equilibrium at an above-zero frequency, i.e. a stable 
coexistence of types occurs. When the interaction coeffi- 
cients of the cooperative strategy are fixed, the equilibrium 
frequency at which it is maintained then depends upon the 
settings of a ss and a cs , i.e. the magnitude of the nega- 
tive density-dependant effects of selfish individuals on them- 
selves and cooperators, respectively. In addition, the within- 
group equilibrium is reached more quickly the greater the 
effects of the selfish type. Although a Lotka-Volterra model 
such as this is typically interpreted at the ecological level as 
representing species interactions, it could also be interpreted 
as a model of allelic competition dynamics within a single 
species group. In particular, coexistence Lotka-Volterra dy- 
namics are analogous to balancing selection for a stable al- 
lelic polymorphism within a group. Conversely, competitive 
exclusion of one species by another is analogous to direc- 
tional selection driving one allele to fixation. The motiva- 
tion for using the language of allelic competition is to fa- 
cilitate comparison with classical group selection models, 
which consider competition between selfish and cooperative 
alleles. 

Our use of the Lotka-Volterra equations in this paper 
should be contrasted fro m their use in community or ecosys- 
tem selection models ( Wilsonl 1992 : Pennl 2003 ). Such 
models do not consider explicit cooperative and selfish types 
in the fashion of traditional group selection models. Instead, 
they examine the complex within-group dynamics that arise 
when a larger number of types are present. These complex 
dynamics can give rise to multiple within-group attractors, 
which can then provide a source of varia t ion in their own 
right upon which selection can act dPennl |2003). By con- 
trast, in this paper we consider simple two-type within-group 
dynamics, where only a single group attractor exists (either 
coexistence or competitive exclusion, as discussed above). 
As far as we are aware, our use of the Lotka-Volterra equa- 



tions to define explicit selfish and cooperative strategies is 
novel. 

Results 

The parameter settings used for the Lotka-Volterra equations 
throughout this paper are shown in Table [T] Changing be- 
tween competitive exclusion and coexistence within-group 
dynamics is achieved by simply switching over the values of 
a ss and a cs , since that determines whether a ss < a cs and 
hence whether competitive exclusion occurs. The values of 
the interaction coefficients in Table Q]produce representative 
within- and between-group dynamics; other settings produce 
the same qualitative trends. In this section, we first present 
results using classical competitive exclusion dynamics, and 
then contrast these to results from the coexistence case. 
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Figure 3: Final group size as a function of time spent repro- 
ducing within groups; initial group size 10 with various % of 
cooperators (competitive exclusion case shown; coexistence 
quantitatively similar). Dotted line shows time at which dif- 
ference in group productivity is greatest. 



Group Selection Dynamics in the Competitive 
Exclusion Case 

The within-group dynamics for a group initialised with unit 
biomass of each type are shown in Figure s |2(a)| and [2(5)1 f° r 
the competitive exclusion case. Initially, both types are in 
their growth phase; their biomass is below the intrinsic type 
carrying capacity of 100. However, the selfish type grows at 
a faster rate, despite the fact that their intrinsic growth rates, 
M, are the same. This is because of the greater negative 
density-dependant effect of the selfish type on cooperators, 
i.e. a cs > a sc . Finally, since a cs > a ss , the coopera- 
tive type is driven to extinction. Furthermore, as Figure [2(b)| 
shows, the proportion of selfish individuals increases mono- 
tonically. Such behaviour is qualitatively identical to that of 
directional within-group selection fo r a selfish allele in clas- 
sical group selection models (e.g. d Wright! 1 19451 : IWilsonL 
1980h ). 
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(a) Biomass of each type. (b) Proportion of selfish type. 

Figure 2: Competitive exclusion within-group dynamics. 

Now let us consider global dynamics under group selec- 
tion in this competitive exclusion case. In order for group se- 
lection to operate through an aggregation and dispersal pro- 
cess, a difference in group size at the dispersal stage must 
exist. Figure [3] illustrates how final group size varies as a 
function of the time spent in the group prior to dispersal, for 
various starting frequencies of cooperators in groups of ini- 
tial size 10. It can be seen from this graph that, using the 
parameters described in Table [T] groups with a greater pro- 
portion of cooperators do indeed grow to a larger size. In 



addition, the results for the coexistence case, where a ss and 
a cs are swapped, are quantitatively similar. These results 
therefore confirm that group selection can in principal oper- 
ate, since there is a variation in group productivity on which 
selection can act. 

To determine the magnitude of the effect of group se- 
lection, the aggregation and dispersal process was executed 
for 5000 iterations, which preliminary experimentation had 
shown to be a sufficient length of time for a global equilib- 
rium to be reached, using the within-group parameter set- 
tings described in Table Q] Equation[3]was iterated 30 times 
in the reproduction stage, while the global population size 
was maintained at 5000. Initial group size was then var- 
ied from 1 to 100 inclusive, while the migrant pool was ini- 
tialised with 50% of each type. The result of this process 
after 5000 aggregation and dispersal cycles is shown in Fig- 
ure |U where 'effect of group selection' on the y-axis is de- 
fined as the difference between the frequency of the selfish 
type at within-group equilibrium and the global frequency 
of the selfish type after 5000 aggregation and dispersal cy- 
cles. Since the within-group equilibrium in the competitive 
exclusion case is the selfish type at 100%, the y-axis equiv- 
ilantely shows the global frequency of cooperators in this 
case. Furthermore, it should also be stressed that the within- 
group equilibrium is the equilibrium that would be reached 
in an unstructured population where there were no groups. 
The y-axis therefore shows the effect that group structure is 
having on the outcome of evolution compared to that in an 
unstructured population. 

There are two points to note from Figure [4] Firstly, in- 
creasing the initial group size decreases the effect of group 
selection, and consequently the global proportion of coop- 
erators. In particular, for small group sizes, the cooperative 
type reaches global fixation (and remains there because we 
do not reintroduce types by mutation). However, for group 
sizes above 10, it is driven extinct. This follows from the 
fact that the between-group variance necessary for group se- 
lection to act is generated by random sampling from the mi- 



Parameter 


Value (competitive exclusion) 


Value (coexistance) 




1.9 


2 


a cs 


2 


1.9 




1 


1 




0.5 


0.5 


K 


100 


100 


M 


0.1 


0.1 



Table 1: Parameter settings of the Lotka-Volterra equation. Note that the only difference between the competitive exclusion 
and coexistence settings is a swapping of the values of a ss and a cs . 
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Figure 4: 'Effect of group selection' (see text) as a function 
of initial group size in the competitive exclusion case. 

grant pool, and therefore rests on the existence of a small 
initial group size, as previously discussed. 

The second, and a key point for this paper, is that the ef- 
fect of group selection rapidly tends to zero as initial group 
size increases. Specifically, above a group size of 10, there 
is no measurable effect at all. Such a result may therefore 
make the idea of group selection acting on randomly formed 
groups seem rather implausible as a significant evolutionary 
pathway. However, the above results only consider the com- 
petitive exclusion case; in the coexistence case, the results 
are somewhat different, as shown in the following section. 

The Efficacy of Group Selection under Coexistence 
Dynamics 

Let us now consider the coexistence dynamics that arise 
from redefining the selfish strategy as a ss = 2 and a cs — 
1.9, i.e. by swapping the interaction coefficients over from 
the competitive exclusion case. Figures [5(a)| and [5 (b)| show 
how the cooperative type is no longer driven extinct at the 
within-group equilibrium. In particular, the change in the 
frequency of the selfish type from an initialisation of 50% 
shows clear balancing selective dynamics resulting in the 
maintance of cooperation at an above-zero frequency. In 
other words, the result is a stable coexistence of cooperative 
and selfish types within a group. 

Group selection dynamics under the aggregation and dis- 
persal process are now as shown in the black curve in Fig- 
ure|6] Crucially, in contrast to the competitive exclusion case 
(shown again in the dotted line), an effect of group selection 
can be seen over the entire range of group sizes examined. 
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(a) Biomass of each type. (b) Proportion of selfish type. 



Figure 5: Coexistence within-group dynamics. 
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Figure 6: Comparing the range of group sizes over which 
an 'effect of group selection' (see text) can be seen between 
coexistence and competitive exclusion dynamics. 



For example, in groups of initial size 50, group selection can 
be seen to still increase the global frequency of cooperation 
above the within-group equilibrium. The significance of this 
observation is that since the within-group equilibrium is the 
same equilibrium that would be reached in an unstructured 
population, these results show that group structure is hav- 
ing an effect on population dynamics across a wide range of 
group sizes. 

Finally, to verify that this result is not an artefact of the 
particular values of a ss and a cs used, the same curves were 
plotted for a variety of other parameters. Figure Q provides 
an example of this, where a ss = 1.99 and a cs = 2 in the 
competitive exclusion case, vice versa for the coexistence 
case. These parameters were chosen since they represent 
stronger within-group selection towards selfish behaviour in 
the coexistence case than in the previous example. Specifi- 
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Figure 7: Demonstrating that the same qualitative trends 
arise where within-group selection towards selfish behaviour 
is stronger in the coexistence case than in Figure [6] Here, 
a cs = 1.99 and a ss = 2 in the coexistence case, vice versa 
for competitive exclusion; all other parameters as in TableQ] 

cally, the within-group equilibrium frequency of the selfish 
type in the coexistence case is 98.04%, compared to 83.3% 
previously. The results in Figure [7] show that while an ef- 
fect of group selection is still seen over a larger range of 
group sizes in the coexistence case, the magnitude of the ef- 
fect is reduced compared to Figure [6] The reason for this 
is that variance in group composition is proportional to the 
frequency of cooperators in the migrant pool in this case, 
and hence to the corresponding within-group equilibrium 
frequency, as discussed in detail in the following section. 

Discussion 

The results in the previous section demonstrate that where a 
stable coexistence of types occurs at within-group equilib- 
rium, an effect of group selection on global frequencies can 
be seen over a much larger range of initial group sizes than 
in the competitive exclusion case. In particular, as group size 
increases in the competitive exclusion case, any measurable 
effect of group selection on the global frequency of cooper- 
ation rapidly tends to zero. By contrast, in the coexistence 
case, some effect on global frequencies is seen over the en- 
tire range of group sizes examined. It must be stressed that 
we do not make a particular claim about the magnitude of 
the effect for large group sizes. Rather, our model implies 
that there is some measurable effect on frequencies over a 
large range of group sizes; how large this effect may be will 
depend on the properties of the natural system under consid- 
eration. However, the fact that any effect of group selection 
still exists over a large range of parameters is significant, 
since it suggests that where within-group dynamics in nature 
are of the coexistence type, some effect of a group popula- 
tion structure may always be acting. 

Coexistence dynamics allow group selection effects to be 
sustained over a larger range of group sizes because of the 
effect of migrant pool frequencies on between-group vari- 



ance. In particular, because group formation constitutes ran- 
dom sampling from the migrant pool, initial between-group 
variance can be approximated by the binomial distribution, 
and is then given by p c p s /S, where p c is the proportion of 
the cooperative type in the migrant pool, p s th e proportion of 
the selfish type, and S the initial group size dWilsonl Il980l 
p. 27). Since p c + p s = 1, it follows that between-group 
variance is proportional to the frequency of the least fre- 
quent type, i.e. variance is maximal when both types are 
of equal frequency, and zero when one type is at fixation. 
Therefore, where one type reaches global fixation then there 
can be no variance and hence no group selection. However, 
in the coexistence case, where one type cannot reach fixa- 
tion, it follows that there must always be some variance and 
hence some possible effect of group selection. The fact that 
the variance is proportional to the frequency of the least fre- 
quent type is illustrated by the difference between Figures [6] 
and |7] where the lower within-group equilibrium frequency 
of cooperators in Figure|7]results in a reduced effect of group 
selection for large group sizes. 

A further observation from Figure |6]is that a gradient to- 
wards an increased effect of group selection also exists over 
a larger range of group sizes in the coexistence case. Specif- 
ically, decreasing group size by a small amount yields an in- 
crease in the effect of group selection for groups of size 20 
in the coexistence case. However, there is no gradient at this 
size in the competitive exclusion case. The significance of 
this is that increasing the effect of group selection increases 
average absolute individual fitness in the population, due to 
an increased global level of cooperation. If group size can be 



partly determined by individual traits (Powers et all 120071) 



then this may provide an adaptive gradient towards smaller 
groups, increased levels of cooperation, and greater fitness. 
In the competitive exclusion case, however, such a gradient 
only exists over a much smaller range of group sizes. While 
numerical experimentation investigating whether either gra- 
dient can be followed by a series of small mutations will 
be the subject of a future study, the results presented here do 
suggest that the concurrent evolution of group size and coop- 
eration is more plausible in cases where a stable coexistence 
of types within groups exists. 

Conclusions 

Any group selection process requires there to be a variation 
in group composition. In aggregation and dispersal style 
models, this variation arises through the random assignment 
of individuals from the migrant pool into groups. Conse- 
quently, it is often suggested that an effect of group selection 
on the global frequency of types will only be seen for very 
small initial group sizes. However, the models on which this 
claim is based typically only consider within-group dynam- 
ics that lead to the competitive exclusion of a cooperative 
type by its selfish counterpart. 

In this paper, we suggest that within-group competitive 



exclusion dynamics may be an unnecessary assumption in 
a number of situations, including the modelling of multi- 
species consortia in bacterial biofilms and during egalitarian 
major transitions. A model has been presented which shows 
that where such coexistence dynamics are present within 
groups, the range of initial group sizes over which an effect 
of group selection can be seen is much larger. Consequently, 
the potential for an adaptive gradient towards smaller groups 
and increased cooperation exists over a much larger range of 
group sizes under coexistence dynamics. This is in sharp 
contrast to the competitive exclusion case, where the ef- 
fects of group selection rapidly reach zero as initial group 
size increases, excluding the possibility of such a gradient 
for a large range of parameters. Our results suggest that, 
where group size can be influenced by individual traits, the 
evolution of smaller groups and increased cooperation is 
more plausible under coexistence dynamics. Such increased 
group cooperation is a vital component of many major tran- 



sitions in evol ution dMavnard Smith and SzathmarvL 1995 
Michodlll999l) . 



We have shown in this paper that the conventional con- 
clusion that group selection effects can only be seen for very 
small groups rests on the assumption that within-group dy- 
namics lead to competitive exclusion. If a within-group co- 
existence of competing types is instead permitted, then the 
range of group sizes over which an effect can be seen is 
much larger. This result follows from the fact that the vari- 
ance in group composition upon which group selection acts 
is dependant not only on group size but also on the frequen- 
cies of types in the migrant pool. In particular, since neither 
type can be driven extinct under coexistence dynamics, there 
will always be some variance in group composition when 
sampling from the migrant pool, which can then be acted on 
by group selection. Thus, since it is not necessary to assume 
that within-group dynamics lead to competitive exclusion, 
this result shows that group selection can operate in a wider 
range of conditions than previously realised. 
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